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Martensitic transformations are first order, diffusion- 
less, displacive, athermal crystal structure transforma- 
tions. They deform the lattice and change the shape of 
the crystal while the relative motion of the atoms is small 
compared to the nearest neighbor distancesi The lat- 
tice deformation of martensitic transformation results in 
a characteristic orientation relationships: Specific vectors 
and planes in the initial lattice are transformed to specific 
vectors and planes in the final lattice. The orientation 
relations can be measured experimentally for many ma- 
terials and provides constraints on the possible atomistic 
pathways of the transformation. Furthermore, marten- 
sitic transformations proceed near the speed of sound, 
complicating a direct observation of the atomic motion.^ 
Despite these conditions on the atomic motion, the atom- 
istic pathway for many martensitic transformations — the 
precise motion of each atom during the transformation — 
remains unknown. 

The problem of finding the most likely pathway for a 
transformation reduces to generating a relevant subset of 
possible pathways and sorting by their energy barriers. 
There are infinitely many unique ways to transform one 
crystal continuously into another. Only a finite subset of 
these, however, results in small strain and small atomic 
motion. Our "pathway generation and sorting" algo- 
rithm produces the set of pathways with small total shape 
change — strain — and small atomic motion — shuffle. The 
pathway generation is related to COMSUBSjS'"' From this 
set, we use the energy barrier of each transformation 
pathway to determine the most likely pathway. 

Titanium's room temperature a phase (hep) trans- 
forms to the high pressure uj phase (3-atom hexagonal 
cell) at approximately 2-9 GPai^ This transformation 
lowers toughness and ductility in Ti alloys. Like all 
martensitic transformations, a^uj occurs at the atomic 
level, so experimentsSiSiSiiSiii can only provide some 
limits on possible mechanisms; moreover, without a sys- 



tematic algorithm for generation, it is impossible to de- 
termine the lowest energy barrier pathway. These prob- 
lems meant that the mechanism for the ol-^uj transfor- 
mation was only recently elucidated by Trinkle el alJ^ 

To apply the pathway generation and sorting algo- 
rithm to Ti a^Lu, we proceed in several steps. We be- 
gin by describing the notation used. We then solve half 
of the pathway generation problem by determining the 
strains for each pathway. Next, we solve the other half 
of the pathway generation problem by determining in- 
ternal relaxations needed for each pathway. Finally, we 
combine the pathway generation methods with the en- 
ergy barrier evaluation methods into a general algorithm 
for producing all relevant pathways and sorting by en- 
ergy barrier, and apply it to the problem of the Ti a— >w 
transformationii^ Finally, the geometry, symmetry, and 
energy barriers of the seven lowest energy barrier path- 
ways are given. 

I. INTRODUCTION 

A. Definition of pathway 

A pathway between two infinite periodic crystals a and 
b is a reversible mapping between the atoms in crys- 
tal a and b. Once each atom in a is identified with an 
atom in b we continuously transform from one crystal 
into another; moreover, the position of each atom during 
this transformation is calculated to determine the energy 
barrier of the transformation. We separate the problem 
of finding mappings between the endpoints — "pathway 
generation" — from following the transformation from a 
to b to determine the energy barrier for finding the low- 
est energy pathway. Once energy barriers are known, we 
choose the lowest energy barrier pathway — the "sorting" 
of the pathways. 
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Pathways for martensitic phase transformations can 
be separated into strain — a global shape change — and 
shuffle — small atomic relaxations. The strain results in 
a macroscopic change of lattice vectors . The shuffles 
produce changes in the atomic basis Xj. 

We restrict ourselves to the mappings that are periodic 
for the entire transformations. This restricts us to homo- 
geneous pathways from infinite crystal a to b. By enforc- 
ing periodicity, we consider only certain transformation 
strains. The shuffles for each pathway are periodic; this 
periodicity allows for a finite search of possible pathways. 
Enforcing homogeneity makes the homogeneous energy 
barrier the relevant "figure of merit" for each pathway. 

We scale the volumes of our crystals so that the vol- 
ume per atom in a is equal to that in b. This simplifies 
the calculation of the strain. For an appreciable volume 
change, a volumetric strain can be added to the resulting 
pathways. 

B. Notation 

A crystal a has lattice vectors written as a matrix [a] 
of column vectors. The volume of the unit cell is det[a]; 
moreover, for the lattice vectors {oi} to be linearly in- 
dependent, det[a] 7^ 0. The atomic basis vectors Xj are 
column vectors in unit cell coordinates; we can represent 
all cartesian points in the infinite lattice using integer 
column vectors fc, R{k;j) — [a]{xj + k). In general, [a] 
and [b] will represent the lattice vectors of crystals a and 
b, respectively; we will also use [a] to denote the lattice. 

A sublattice [A] of a lattice [a] is a lattice where every 
point in [A] is also in [a]. The lattice vectors [A] are 
called the supercell. In order for [A] to be a sublattice of 
[a], there must exist an non-singular integer matrix [n] 
such that [A] = [a\[n]. This condition guarantees that 
any point in [A] is also a point in [a]. The size of the 
supercell [A] is the ratio of the number of points in [a] to 

[A] ; it is det[n]. Generally, we will write [A] and [B] for 
a supercell of [a] and [6], respectively; the integer matrix 
[n] will relate [a] to [A] , and [m] will relate [b] to [B] . 

Atoms in crystal a are mapped into new positions in 
the crystal A. The basis for A is {u°}; each atom in A is 
identified with some atom in a: [A]uf — [a] (xj+k) , where 
fc is a column vector of integers. This means that uf = 
[n]~^{xj + fc) for some integer vector k. We require that 
each component of u° be in [0,1), so there are exactly 
det[n] possible vectors fc. This means that if a has Na 
atoms, then A will have Na det[n] atoms. 

A pathway is described as a supercell pair [A] and [B] 
with TV atom positions Ui and shuffles Si. Each super- 
cell relates to its crystal by — [a]ixj + k) and 

[B] {ui -\-5i) = -l-fc') -I- (constant shift), where fc and 
fc' are integer vectors and the constant shift is a uniform 
translation of every point in the lattice. These equations 
connect the atom positions Ui and shuffles 5i to the un- 
derlying crystals a and b. Moreover, the supercells [A\ 
and [B] are related by a symmetric strain tensor e and a 



rotation matrix <I>, 

e[A]=<^[B]. (1) 

The strain e changes the lattice vectors \Al\ into those of 
[B] after a rotation. This corresponds to a change in the 
lattice vectors [a] since [A\ = \a][n]. Moreover, the strain 
is determined entirely by the supercells. 

These equations allow us to write each pathway as 
P{[n],[m];ui,5i), where [A\ = [a][n] and [B] = \b][m]. 
The advantage of working in terms of the integer ma- 
trices [n] and [m] is the convenience in enumerating all 
possible supercells of interest. We will often talk about a 
supercell pair [n] and [m] and write P{[n]^ [m]; *) to indi- 
cate that atom positions and shuffles have not yet been 
determined for this pathway. 

We determine the strain, orientation relations, dis- 
tance that atoms move, distance of closest ap- 
proach to other atoms, and common subgroup from 
P{\n],[m]]Ui,5i). The strain e is determined purely 
by the supercell pair P([n], [m]; *). The orientation 
relations — which vectors and planes in the lattice [a] are 
transformed into vectors and planes in the lattice \b] — 
are also determined from P{[n\, [to];*). The shuffle in- 
formation is sufficient to determine a continuous linear 
simultaneous transformation from a to b; we use this 
transformation to calculate the distance that each atom 
moves and the distance of closest approach. 

We determine whether two pathway representations P 
and P' are different by comparing their supercell pairs 
and their shuffles. At each step, equivalent pathway rep- 
resentations are removed. The details of these compar- 
isons are in Section Hi CI for the supercell pairs and Sec- 
tion ^^D] for the shuffles. 



II. SOLVING THE STRAIN PROBLEM 

The first step in enumerating possible pathways is to 
find the set of supercell pairs P{[n]^ [to]; *). We find all 
supercell pairs of a given size N where the strain e is lim- 
ited by a cutoff £max- We begin by enumerating possible 
unique sublattices of a given size, and reduce this list by 
the lattice symmetry. We then revisit Eq. and solve 
for the strain e and rotation $. We restrict strains to a 
cutoff Emax to producc a finite search algorithm enumer- 
ating all possible supercell pairs. Some aspects of this 
problem were originally considered in Lomer's calcula- 
tion of the orientation relations in U fi^aw^ 



A. Unique sublattices 

We say two lattices x and y with lattice vectors [x] and 
[y], respectively, are equivalent if every point in [x] is a 
point in [y] and vice versa. This mirrors the definition of 
sublattice from Section [IBI i.e., [x\ and [y] are equivalent 
if and only if [x\ is a sublattice of [y] and [y] is a sublattice 
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of [x]. This means there are two non-singular integer 
matrices [i] and [j] such that [x] — [y] [i] and [y] — [x] [j] . 
Combining these equations shows that [i] and [j] must be 
unimodular, that is, det[i] = ±1 and det[j] — ±1. This 
gives another condition for equivalence: There exists an 
integer unimodular matrix [/] such that [x] — [y] [I] . 

We can apply this definition to supercells as well; two 
supercells, given by [n] and [n'] are equivalent if and only 
if there exists an integer unimodular matrix [I] such that 
[n] = This also means that equivalent supercells 

must be of the same size; i.e., |det[n]| — |det[n']|. We 
use unimodular matrices to iterate over possible super- 
cell representations for a sublattice; this is equivalent to 
forming alternate supercells that have the same size. 

From each set of equivalent supercells, we pick one 
representative called [n]; all other members of the set 
can be generated using integer unimodular matrices [I] 
as [n] [I] . The representative sublattices are given by an 
upper-triangular matrix [n], where n,j > for all 

fiij < fiii for j > i, and Yl^^i nu = N, where N is the 
size of the sublattice det [n] . It is straightforward to show 
that for a given N there are a finite set of supercells in 
this form. If two matrices [n] and [fi'] have the above 
form, they are equivalent if and only if they are equal 
(see Appendix D of 14]). 

Using the form for representative supercells we enu- 
merate all possible unique sublattices for [a] and [b] . The 
set of unique sublattices for [a] is {[n]} and for [b] is 
{[m]}. We restrict ourselves to pathways of a given size 
N; this requires that iVa det [n] ^ N ^ A^f,det[m], where 
Na and are the number of atoms in crystal a and b, 
respectively. 

To further reduce the set of representative supercell, 
consider the set of all symmetry operations Go on [a] 
where at least one lattice point is mapped to itself; this 
is the point group of a. We can represent each member of 
the point group Ga with a matrix ga that operates on the 
cartesian coordinates of [a]. In order to be a symmetry 
element, ga must map the lattice [a] back onto itself: 
There exists an integer unimodular matrix ga such that 
(7a [a] = [a\ga- Since the lattice a is left invariant under 
the operation of ga, two initially different sublattices A 
and A' may be equivalent by ga- 

To determine the set of equivalent sublattices of size 
N for a given lattice a, we start with the set of unique 
representations {[n]} and reduce it by symmetry. We say 
that two of our (initially unique) representations [n] and 
[h'] are equivalent if there exists a symmetry element ga 
(with integer unimodular matrix ga) and some integer 
unimodular matrix [i] such that ga [n] = [n'] = [n'] [i] . 
Thus, our initial list of unique sublattices of size N may 
be further reduced by the symmetry of lattice a. Such a 
reduction can also be applied to the sublattices of lattice 
b. 



B. Calculating and limiting the strain 

Given two sets of sublattice representatives {[n]} and 
{[m]} for a and b, respectively, we find all strains that 
transform from one sublattice into the other. We limit 
the allowed strain by a maximum cutoff £max, which pro- 
duces a finite list of possible supercell pairs. We begin by 
solving the strain Eq. |^ for the general case. We then 
limit the allowed strain, and translate that limit into a 
subset of allowed supercell pairs. 

We take a supercell representative [n] for a and [fh] 
for b and substitute into Eq. to get an equation for 
the possible strains from [h] to [fh]. The supercells are 
[A] — [a] [h] [i] for some integer unimodular matrix [i] , 
and [B] — [b][fh][j] for some integer unimodular matrix 
[j]. The equation for e and $ is £[a][fi][i] = $ [6] [to] [j] , 
which can be simplified by right multiplying by [i]^^, 

e[a][n]=m[rn][lh (2) 

where [/] = [j][i]~"'^ is an integer unimodular matrix. 

To find all possible strains, we will solve Eq. ||2J) for all 
integer unimodular matrices [I] to find the supercell pairs 
P([n], [rh][l];*). We right multiply Eq. @ by [fi]-i[a]-^ 
to get £ — ^[b][fh][l][h]~^[a]~^ . We define the matrix 
C{[1]) = [6][TO][;][n]-i[a]-i; then e = ^C{[1]). The strain 
tensor e is symmetric, though C{[1]) and $ are not. We 
left multiply each side of the equation by its transpose to 
get 

where the second equality is true because e is symmetric 
and $ is a rotation matrix, so = Thus, e is 

the square root of C([Z])'^C([Z]), a symmetric matrix. We 
calculate the square root by diagonahzing C([Z])'^C([Z]) 
and writing it as Q"^cQ where is a rotation matrix 
and c is a diagonal matrix of the eigenvalues. Then e = 
e'^y/cQ, and ^ are the diagonal strain elements of e; 
because there is no volumetric change, the product of 
the diagonal strain elements are unity. Once e is known, 
<i> = e[a][n][/]"^[TO]~^[6]~^ This method may also be 
considered an extension of the method of magic strainsi^^ 
Though every integer unimodular matrix [/] has an as- 
sociated strain tensor e, most of these will be large. We 
define a strain limitation in terms of the cutoff Emax > 1; 
a strain e is within our cutoff if for all non-zero vectors 

£mLl«l <\£V\< Emaxlw]- (4) 

This condition is equivalent to requiring that all the 
eigenvalues of our strain matrix are between e~ax s-'^d 
£max- We consider only those [I] that will produce a strain 
within our £max cutoff. 

Eq. 1^ relates each supercell lattice vector Ai to a 
supercell lattice vector Bf, this translates Eq. Q to a 
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condition on each Bi. We substitute Ai for v in Eq. Q 
noting that eAi is Bi, so that for each i, 

e^Ll^.l < l^.l <e„,ax|i',|. (5) 

Thus, if {Bi} doesn't satisfy Eq. (jSJ, then e will never 
satisfy Eq. Hence, only Bi which lie in the annulus 
between ej^J^^l^il and £max|^i| need to be considered; 
this is a finite list for each Bi. 

We translate Eq. lO into a condition on the allowed 
unimodular integer matrices [/]. For each i, we con- 
struct the (finite) set of integer vectors {A^*''} such that 

Bi = [b][m\xf satisfies Eq. (0. This set will be different 
for each i, as well as for each supercell [n] and [m\. To 
construct all of the possible unimodular matrices [I] , we 

check to see if the matrix j2, Js)] = (a^^^ A^-^^ A^^^^ 
has determinant f. If it does, we determine e for the 
supercells [n] and ['rh][l{ji, j2, js)], and make sure that 
e is within the cutoff emax- To efficiently construct the 
{Aj*'} we choose representations [n] and [rh] that have 
the shortest possible lengths in cartesian coordinates for 
unit cells [a] and [b], respectively. This makes both \Ai\ 
as short as possible for each [n] , and makes the set of A^*' 
to check as small as possible. 



As an example, consider the cubic crystal class; it has 
only three unique elastic constants Cn, C12, and C44. 
We write the energy for a general strain a as: 

Uaie) = ^Ciiiej +el + ej) + Ci2(eie2 + 6263 + egCi) 

+ lCii{el + el + el). 

To determine if two strains are truly identical, we check 
that the three combinations + e| + e§, 6162 + 6263 + 
6361, and 64 + e| + e| are the same for e and e' . Similar 
expressions can be derived for the other crystal classes. 

III. DETERMINING SHUFFLES 

Each supercell pair [m]; *) represents the strain 

of the transformation; to complete the pathway descrip- 
tion requires knowledge of the local atomic motion: the 
atom positions Ui and shuffles 6i. We find the shuffles 
using a systematic approach to enumerate all possible 
mappings between the atoms in the supercell [A] and 
the atoms in [B]. We use geometric restrictions to re- 
duce the set of possible pathways to those that require 
small atomic motions and do not have artificially small 
approach distances during the transformation. 



C. Uniqueness of supercell pairings 

Two different supercell pairs P([n], [m][Z]; *) and 
P([n'], [TO'][r]; *) are equivalent if they produce the same 
mapping from a to b. This first requires that [n] = [h'] 
and [fh] = [m'] since supercell representatives are equiv- 
alent if and only if they are equal. However, [I] and [I'] 
need not be equal to give equivalent mappings if they 
produce equivalent strains from a to b and vice versa. 

Two strains e and e' on a crystal a are equivalent if for 
arbitrary elastic constants C° the elastic energies Ua{s) 
and Uaie'), 

1 ^ 

are equal. The strains e,; are related to the strain matrix 
e by 

(1 + ei iee ^65 \ 
ieg 1-1-62 564 . 
265 ^64 1 + 63 / 

The symmetry of the elastic constants C°j is given by the 
classic of a; there are 11 unique three-dimensional crystal 
classes, and 5 unique two-dimensional crystal classes. We 
perform the same test for the strains on crystal b. This 
test takes into account the point group symmetry of the 
lattices a and b. 



A. Populating cells 

For a given supercell pair P([rt], [m]; we popu- 
late the [A] supercell with the atoms from a and the 
[B] supercell with the atoms from b using the formula 
from Section II BI The atom positions in [A] are 
uf = [n]^"'^(a;^ + k) for some integer vector fc; similarly, 
Ui — [m]~^{x'j + k'). We restrict all of the components 
of Ui and to be in [0, 1), so there are exactly N atoms 
from each crystal. 

The shuffle vectors 6i connect uf to the atoms u'j in 
the periodic supercell as the strain e connects [A] to [B] . 
Because of the periodicity of the supercells, the mapping 
can connect an atom i to an periodic image of an atom 
j that is not inside the original supercell. Suppose we 
have our mapping, and each uf is transformed into some 
^i(i) where hj is an integer vector; our shuffle vectors 
Si are 6i = u^^.^j -|- hj - uf. 

Our problem then is to find the mapping for each atom 
i in [A] to an atom j in [B]. Moreover, we look for map- 
pings that result in small atomic motions and do not 
bring two atoms close together, both of which produce 
large energy barriers. 

B. Center of mass 

We require that the shuffles do not change the center 
of mass of the crystal during the transformation: Si = 
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0. This is enforced by shifting all of the atoms by a 
constant vector ^com- Then, Si = u^j-^j +ScoM + hj—uf. 
We solve for (5com by summing Si over all N atoms 

^Si=0 = + '^COM + hj - u"), 

i i 

i i 

Then 

5coM = ^E«-"')-^E^.- 

i j 

The only term in this equation that depends on the spe- 
cific mapping is the sum over integer vectors hj. Thus, 
every center of mass shift (5com has the form 

ScoM = ^ - ^4) + ^(fci, fc2, fca), (6) 

i 

where ki = Q . . . N — 1. We need not consider shifts 
larger than this, as shifts by integer vectors only translate 
all of the atoms u\ from one periodic supercell image to 
another. 

Because every mapping will have a center of mass shift 
given by Eq. ©, we systematically attempt each shift 
ScoM on the atoms m^, and make the mappings from 
uf to keeping only those where ^ ■ Si = 0. We still 
consider mappings where uf moves to + hj , but restrict 
the hj components to all be either -1, 0, or 1. This 
prevents an atom from moving across an entire supercell. 

The center of mass shifts i5com are reduced by con- 
sidering only shifts that stay within the unit cell of a and 
b. There are some shifts k/N where [m\k/N is all inte- 
ger; this corresponds to shifting all of the atoms in b by 
a lattice vector of [6]. This will create no unique path- 
ways, but rather previously considered pathways with a 
permutation of the atom indices j . Similarly, some shifts 
correspond to shifting all of the atom by a lattice vector 
of [a]; these need not be considered either. We reduce 
the set of possible (5com to those that are within a sin- 
gle unit cell of [a] (defined by [n]) and [b] (defined by 
[m]). Depending on the orientation of [n] and [m], this 
reduces the number of center of mass shifts to somewhere 
between N^/{NaNb) and / ma.x{Na, Nb). 



C. Mapping atoms for a given center of mass shift 

For each supercell pair P([n], [m]; we loop over the 
set of possible center of mass shifts ^com and enumerate 
all possible mappings for each center of mass shift. This 
becomes a combinatoric problem; to reduce the exponen- 
tial scaling, we examine only mappings where no atom 
moves "far," given by a shuffle cutoff ^max- We further 
reduce the set of possible mappings to those where atoms 



do not approach closely and where the total shuffle for 
all atoms is small. 

We define a metric d{x — y) between two vectors x and 
y in supercell coordinates 



d{x -y)^ \l2\^A]{x - y)P + -\[B]{x - (7) 

This distance function is symmetric in the supercells [A\ 
and [B], is zero if and only if x = y, and obeys the 
triangle inequality. We use this function to compare the 
atom positions uf and even though they are defined 
from different supercells; only pathways where all atoms 
move less than (Smax are considered. 

Given our distance function, we calculate the distance 
d{u'j + h — uf) for all i and j and shifts h. We construct 
the possible shifts h for each i and j on a component- 
by-component basis to minimize the number of h^s to 
consider. If the l^^ component of — uf is less than 0, 
we use /i; = or /i/ = 1; if instead it is greater than 0, 
we allow ft,/ = or ft,/ = —1. In this way, we construct 2^ 
shifts h for each pair This produces a table of N by 
8N entries of distances. 

For each atom i we construct the set of atoms C{i) = 
{u'j + ft} where each atom is within Smax of uf. These 
are the allowed atom identifications for each i. Our com- 
binatorial problem is on the sets C(z); no two atoms i 
and i' may map to the same atom j. The set of possible 
pathways is found by iteration, performed using recur- 
sion. For i — 1, we pick in turn one element of C(l); this 
will be j(l). For each i > 1, we remove any entries of j(l) 
from C{i); call the new sets C(i)(i). We then pick in turn 
one element of C(i)(2); this is j{2). For each i > 2, we 
remove any entries of j(2) from C(i)(i); call the new sets 
C(i,2) (*)• We repeat for successive i until i = N, at which 
point we have an entire mapping If at any point a 
set C(i_..._fc)(i) for i > k becomes empty or we exhaust 
all the possible entries in a set we go back and make a 
different choice for some smaller i. This will produce all 
possible mappings consistent with the initial sets C{i). 

After each possible pathway is constructed, we cal- 
culate the Si and only keep pathways where Si is 
zero. Because we consider all possible Scom in our "outer 
loop," we are guaranteed to enumerate all possible path- 
ways within our i5niax limit. We store the possible path- 
ways for each center of mass shift, and then check that 
each pathway is unique. 

D. Uniqueness of pathways 

Because each atom is indistinguishable, we define 
equivalence of pathways in terms of the local environ- 
ment that atoms will see along the pathway. We believe 
our definition for equivalence is sufficient, but have not 
proven such. This was tested by comparing against the 
pathways produced in the common subgroup method^ 

We use three tests to determine if two shuffle sets Si 
and S'^ are equivalent. First, the sets of shuffle magnitudes 
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must be equal: 



{md,i\[B]5,\),. 

{{\[A]S[l\[B]S[\),. 



.A\[A]SMA[B]Sn\)} = 
■,(l[^]^^vU[*l)}- 



If this is true, then we next require that distance of closest 
approach during the transformation — dca — match. We 
calculate dca by simultaneously linearly interpolating all 
atoms from their initial to final positions using a single 
variable x while straining the cell by 1 + a;(£ — 1). If 
the distance of closest approach is the same, then we 
check the nearest neighbor distance for each atom in the 
half-strained, half-shuffled "intermediate" supercell. If 
those sets are the same, then we say that we have two 
equivalent pathways. 

These tests only comprise a set of necessary conditions, 
not sufficient, for equivalence; this should not be seen as 
a severe limitation. First, two inequivalent paths marked 
as "equivalent" would be similar insofar as the trans- 
formation would produce similar local environments for 
the atoms. It's not entirely clear there would be a very 
large difference in energy as the crystal transformed us- 
ing either pathway. Second, our choice of tests has been 
checked by comparing to an alternate method of atom 
identification using the Wyckoff positions in the strained 
"intermediate" lattice; both methods produced identi- 
cal lists of unique pathways for many different supercell 
choices. While we don't know that the Wyckoff tests are 
not also limited, they should not have the same limita- 
tions, and so it is very likely that the limitations of each 
must be very small. Finally, the real test of equivalence is 
that two pathways are equivalent if for each atom in the 
first pathway, there is one and only one atom in the sec- 
ond pathway that has the same local environment during 
the transformation. Such a test would require a compli- 
cated A^! search, and would ultimately make the test for 
uniqueness computationally intractable. We believe our 
set of conditions is a useful subset of the full equivalence 
test. 



E. Reduction by total shuffle magnitude and 
distance of closest approach 

We remove pathways with small distances of closest 
approach and large root mean square (RMS) shuffle mag- 
nitudes because their energy barriers will be high. The 
majority of pathways generated for a given Jmax are en- 
ergetically unlikely because they require two atoms to 
come very close together, or the majority of the atoms 
to move a large distance. For each supercell pair, we re- 
tain a subset of pathways by examining their RMS shuffle 
magnitudes 



<5rms 



1 ^ 



\ 



1 ^ 



and the distance of closest approach for each. We 
find the largest distance of closest approach for the 



supercell pair P{[n],[m]; *), and set the minimum al- 
lowed distance of closest approach to be min(dca) = 
max((ica) — 0.1 A. For the supercell pair, the maxi- 
mum allowed RMS shuffle magnitude is max((5RMs) = 
max (V2min(^RMs), feMS (best dca pathway)). 

These two rules (a) reject many of the poor candidates 
for each supercell pair, (b) while ensuring that at least 
one pathway for each supercell pair is examined, even if 
it has a small distance of closest approach. If a supercell 
pair has a "good" solution (small Sbms and large dca.) we 
check only a few other possible pathways; if the supercell 
pair does not have any "good" solutions, we check many 
possible pathways for that pair. 



IV. CALCULATING ENERGY BARRIERS 

Once a set of pathways P{[n\, [rn\[l];ui,Si) is known, 
we determine which pathways are the most probable for 
a given material. The sorting criterion we use is that 
the pathway with the lowest energy barrier should be 
the most likely to occur during a homogeneous transfor- 
mation. Most martensites are believed to transform by 
heterogeneous nucleation, and there are some caveats in- 
volved in translating the homogeneous pathway into a 
heterogeneous pathway. However, those limitations are 
material specific, and so cannot be solved for the general 
case; these issues have been addressed for titaniumjiSi 

Calculation of the energy barrier of a pathway requires 
(a) an atomic interaction potential and (b) a method for 
determining the energy barrier of a transformation. The 
latter can be solved accurately with some computational 
effort, or approximately with less effort. Because ulti- 
mately only the pathways with the lowest energy barriers 
are important and not the many higher energy pathways, 
we use approximate barrier calculations for the vast ma- 
jority of the pathways. This leaves the accurate calcula- 
tion of the barrier to a subset of the most likely pathways, 
reducing computational effort. 

We use three methods to calculate the energy bar- 
rier, from most accurate and computationally inten- 
sive to least. The first method — the nudged elastic 
band method^^ — gives the accurate transformation bar- 
rier. The landscape barrier uses a reduced phase space of 
one strain and one concerted shufffe variable. The lack of 
relaxation makes the landscape barrier an overestimate 
of the nudged elastic band barrier. The elastic barrier is 
an approximation to the landscape barrier, constructed 
only from the strain. The lack of atomic "stiffness" from 
phonons makes the elastic barrier an underestimate of 
the landscape barrier. Details of the landscape and elas- 
tic barriers have been published elsewhere^ 



V. TITANIUM a^uj RESULTS 

The pathway generation and sorting algorithm was ap- 
plied to the Ti a— >a; transformation, and a new pathway 
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emerged with an energy barrier lower than all others by 
over a factor of four."- Herein we provide a detailed study 
of the methodology used and the remaining low energy 
barrier pathways studied. 



A. Computational details 

Parameters for the pathway generation algorithm are 
chosen to ensure sampling of all relevant pathways, and 
the total energies are computed to yield accurate energy 
barriers. The 2-atom unit cell of hep and 3-atom unit 
cell of omega require that all supercell pairs contain a 
multiple of 6 atoms; so we use = 6 and 12. We re- 
quire om diagonal strains be less than Emax = 1.333; 
for comparison, the Burgers pathway for hep to bcc re- 
quires a diagonal strain component of only l.li^^ We 
keep only pathways with elastic barriers less than Ecut = 
100 meV/atom. From those supercell pairs, we construct 
pathways with a maximum shuffle magnitude (5inax less 
than 2.0 A; this is comparable to the interatomic dis- 
tance in a and w, and, as seen later, is much larger than 
the shuffles in our lowest energy pathway. Each super- 
cell pair produces multiple pathways with shuffles smaller 
than this value. As discussed in Section Fill El we keep 
pathways with closest approach distance dca within 0.1 A 
of the largest dca- We reject any pathways with an RMS 
shuffle (5rms which is larger than both the Sjxms of the 
best dca pathway found and V2 of the smallest (5rms- 
This gives us a set of pathways with good distance of 
closest approach, RMS shuffle, and guarantees at least 
one pathway for each supercell pair. 

We calculate total energies for the landscape energy 
barriers using a tight-binding (TB) model. The TB cal- 
culations are performed using the molecular dynamics 
code OHMMsiSi and use Mehl and Papaconstantopoulos's 
functional form22i2i with parameters refit to reproduce 
full-potential density-functional total energies for hep, 
bcc, fee, omega, a.nd sc to within 0.5 meV/atom (see 
Appendix B of [iJI). For each mechanism, we use a k- 
point set equivalent toal2xl2x9orl2xl2x8 grid in 
the original hep lattice. This k-point mesh is converged 
to within 0.1 meV/atom against a 18 x 18 x 12 k-point 
mesh for the three lowest energy barrier pathways. A 
Fermi broadening of 63 meV (5 mRyd) was used with 
this k-point mesh to ensure a smooth electronic density 
of states. 

We calculate total energies and forces for the nudged 
elastic band barrier using carefully converged ab ini- 
tio calculations, performed with VASP.^^'^^ Vasp is a 
plane-wave based code using ultra-soft Vanderbilt type 
pseudopotentials^i as supplied by Kresse and Hafncr.^^ 
The calculations were performed using the generalized 
gradient approximation of Perdew and Wangi2& We in- 
clude 3p electrons in the valence band and use a plane- 
wave kinetic-energy cutoff of 400 eV and a 7 x 7 x 7 
k-point mesh to ensure energy convergence to within 
1 meV/atom. We relax the atomic positions and the 



6 atom 12 atom 

Sorting supercell supercell 
step: pairs pathways pairs pathways 



1. 


£max 


6 




128 




2. 




4 




56 




3. 




4 


55 


56 


3790 


4. 


dca, SkMS 


4 


6 


56 


971 


5. 


-^landscape 


2 


2 


3 


5 



TABLE I: Number of supercell pairs, pathways generated, and 
cutoff parameters in Ti hep to omega pathway search. Each 
consecutive step builds on the selection made in the previous 
step: 1. We generate an initial set of 134 supercell pairs with a 
given Emax cutoff of 1.333. 2. We choose the 60 supercell pairs 
that have an elastic barrier less than Ecut of 100 meV/atom. 
3. The possible shuffles for each supercell pair within 5max of 
2 A are determined to give 3,845 possible pathways. 4. For 
each supercell pair, we define a minimum allowed distance 
of closest approach (dca), and a maximum root mean square 
(RMS) shuffle magnitude ((5rms), and keep the 977 pathways 
within those limits. 5. Finally, the 7 pathways with tight- 
binding (TB) landscape barriers less than 90 meV/atom are 
kept, and their true energy barriers calculated. 



unit cell shape and volume until the atomic forces are 
less than 20 meV/A and the stresses are smaller than 
20 MPa. 



B. Generated pathways 

Table m summarizes the resulting number of supercell 
pairs and pathways generated and the cutoffs made at 
each step in the algorithm. The nudged elastic band 
energy barrier for the final seven pathways is calculated, 
and the pathway with the lowest barrier is TAO-1.^^ 

We generate an initial set of 134 supercell pairs with 
diagonal strain components less than Emax and reduce 
it to 60 pairs with elastic barriers less than Ecut of 
100 meV/atom. The three-dimensional surface of possi- 
ble diagonal strain values Si satisfying e~ax ^ ^ £max 
and £i£2£3 = 1 has sharp corners at the boundary; be- 
cause of this, we generate more supercell pairs than are 
used to ensure that we produce all supercell pairs less 
than our elastic barrier cutoff. This procedure is analo- 
gous to inscribing a circle inside of a regular hexagon. 

Using larger supercells fails to produce any supercell 
pairs with strains smaller than the smallest strains al- 
ready found. The possibility for better pathways for 
larger supercell sizes is unlikely as there are no 18- or 24- 
atom supercell pairs found with strains smaller than the 
smallest strain found using 1 2-atom supercells. We keep 
the 12-atom supercell pairs because there is a 12-atom 
supercell pair that has smaller strains than the small- 
est strain for a 6-atom pathway; three of the seven low 
energy pathways have this smallest strain valuei-^ 

The 60 supercell pairs generate 55 6-atom pathways 
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and 3790 12-atoni pathways with ^max of 2.0 A; this set 
is reduced to 6 6-atoiii and 971 12-atom pathways using 
the closest approach distances and RMS shuffle magni- 
tudes as described in Section IIII El The reduction for 
the 6-atom pathways is more severe because each of the 
supercell pairs has a "good" solution (small RMS shuf- 
fle magnitude and large closest approach distance) so we 
examine fewer alternative pathways. Many of the 12- 
atom supercell pairs do not have "good" solutions, so we 
examine more alternative pathways. A check of this re- 
duction for good solutions was performed by calculating 
the landscape barrier of the "best" excluded pathway for 
the TAO-1 6-atom supercell pair; it has a TB landscape 
barrier of 160 meV/atom, compared to 43 meV/atom for 
TAO-1. 

The 977 pathways are then reduced to seven: three 
with landscape barriers less than 80 meV, and four with 
landscape barriers between 80 meV and 90 meV. The 
three lowest pathways are TAO-1 (6-atom), TAO-2 (6- 
atom), and Silcock (12-atom); the latter is a pathway 
already proposed in the literature^ The remaining six 
pathways are new; though TAO-1 and TAO-2 are related 
to Usikov and Zilberstein's a—ff3—fUJ pathway through 
an intermediate bcc structured We calculate the nudged 
elastic band barrier for Silcock and TAO-1 through TAO- 
6. 

Tablenigives details for the seven lowest-barrier path- 
ways of interest. Along each pathway, the symmetry is 
a subgroup of both crystals a and b. Using COMSUBSj^ 
we obtained the maximal symmetry of this subgroup for 
each pathway. In the table, we list this maximal com- 
mon subgroup along with atomic positions in the setting 
of this subgroup. This gives a complete description of 
the mapping of atoms in a onto atoms in b. For exam- 
ple, in the TAO-1 pathway, the symmetry of the crystal 
is monoclinic (space group 15 C2/C) as it transforms 
from a to to. The atoms are at the Wyckoff (e) and (/) 
sites of that space group. During the transformation the 
atoms at the (e) site {0,y, 1/4) move from (0, 7/12, 1/4) 
to (0, 1/2, 1/4) and the atoms at the (/) site {x, y, z) move 
from (1/3,5/12,1/2) to (1/3,1/2,1/6). 

It is important to note that in the nudged-elastic band 
method, the atomic positions are relaxed under triclinic 
symmetry: no point symmetries are assumed or enforced. 
We find the resulting symmetry along each pathway to be 
nonetheless equal to the maximal symmetry given in Ta- 
blem In each of these cases, we conclude that the lowest 
barriers occur for pathways with maximal symmetry. 

Table Hill summarizes the energy barriers for the seven 
lowest pathways of interest. The TB NEB barrier is 
bounded above by the landscape barrier and below by 
the elastic barrier. The difference between the TB NEB 
barrier and the ab initio NEB barrier is due primarily to 
the latter's cell size and shape relaxations. 

Figure ^ shows the energy barrier for the seven lowest 
energy barrier pathways. The plots are grouped accord- 
ing to the supercell pairs in each pathway. The difference 
in energy barrier among pathways with the same super- 
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FIG. 1: Energy barrier at zero pressure for the seven low- 
est energy pathways, using ab initio NEB at 16 intermediate 
states. The pathways are grouped by their strains; each set 
have related supercell pairs. The effect of different shuffles on 
the energy barrier can be seen by comparing TAO-2 to TAO- 
5, and Silcock, TAO-3, and TAO-4 to each other. The latter 
three pathways have the smallest amount of strain, but the 
homogeneous barrier is smallest for TAO-1. 



cell shows the importance of theory in understanding 
the microscopic mechanism — measurements of the final 
Lu orientation relative to the a can only rule out possible 
pathways. 

Figure |21 illustrates a simple geometric picture of the 
low barrier TAO-1 pathway. The basal plane of o; is a 
series of hexagons surrounding individual atoms; these 
hexagons then break into two 3-atom pieces. Each 3- 
atom piece swings in opposite directions out of the basal 
plane, and connects with a mated 3-atom piece above or 
below. These form the honeycomb sublattice of uj; the 
remaining "unmoved" atoms form the A sublattice in uj. 
This transformation transforms the (0001) plane of a into 
the (0111) plane of oj. 

Of the remaining six pathways, TAO-2 and Silcock are 
both related to existing pathways in the literature, while 
TAO-3 through TAO-6 are related to the other three 
pathways. TAO-2 has the same endpoints as Usikov and 
Zilberstein's variant II pathway^ but avoids the inter- 
mediate bcc phase. The Silcock pathway is identical to 
the original published pathway. Both TAO-3 and TAO- 
4 use the Silcock supercell pairing, though with different 
shuffles. TAO-5 is a ceU doubling of the TAO-2 superceU 
pairing with different shuffles. Finally, TAO-6 uses the 
same a supercell as Silcock, but connects it to a different 
LU supercell. 
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Wyckoff positions 


Path 


Subgroup 




Supercell lattice 


m 


label 


a 




TAO-1 


15 C2/c 


a 


(-1,-1,2), (1,-1,0), (1,1,1) 


2 


(e) 


0, 7/12,1/4 


0,1/2,1/4 


(6-atom) 




U! 


(-2, -1,2), (0,-1,0), (0,0,2) 


4 


(/) 


1/3, 5/12,1/12 


1/3,1/2,1/6 


TAO-2 


11 P2i/m 


a 


(1,0,0), (0,0, 1),(0, -3,0) 


2 


(e) 


1/3,1/4,7/9 


1/2,1/4,3/4 


(6-atom) 






(0,0,-1), (-1,0,0), (1,2, -2) 


2 


(e) 


1/3,1/4,1/9 


1/3,1/4,1/12 










2 


(e) 


1/3,1/4,4/9 


2/3,1/4, 5/12 


Silcock 


57 Pbcm 


a 


(-1,0,0), (-3, -6,0), (0,0,1) 


4 


(d) 


3/4, 23/36, 1/4 


1/2, 5/8, 1/4 


(12-atom) 






(0,0, 1), (2, 4,0), (-1,0,0) 


4 


(d) 


3/4, 11/36, 1/4 


1,7/24, 1/4 










4 


(d) 


3/4, 35/36,1/4 


1,23/24,1/4 


TAO-3 


6 Pm 


a 


(1,0,0),(0,0,1),(0, -6,0) 


1 


(a) 


1/3,0,1/18 


0,0,0 


(12-atom) 






(0,0,-1), (-1,0,0), (2, 4, -3) 


1 


(a) 


1/3,0,2/9 


0,0, 1/6 










1 


(a) 


1/3,0,7/18 


1/2,0,1/3 










1 


(a) 


1/3,0,5/9 


1/2,0,1/2 










1 


(a) 


1/3,0, 13/18 


1/2,0,2/3 










1 


(o) 


1/3,0, 8/9 


0,0,5/6 










1 


(6) 


2/3, 1/2, 1/9 


1/4, 1/2, 1/12 










1 


(&) 


2/3, 1/2,5/18 


1/4, 1/2, 1/4 










1 


(&) 


2/3,1/2,4/9 


5/4, 1/2,5/12 










1 


(6) 


2/3,1/2,11/18 


3/4, 1/2, 7/12 










1 


(b) 


2/3,1/2,7/9 


3/4, 1/2, 3/4 










1 


{b) 


2/3,1/2,17/18 


3/4,1/2,11/12 


TAO-4 


11 P2i/m 


a 


(1,0,0),(0,0,1),(0, -6,0) 


2 


(e) 


1/3, 1/4, 13/18 


-1/8, 1/4, 17/24 


(12-atom) 






(0,0,-1), (-1,0,0), (2, 4, -3) 


2 


(e) 


1/3,1/4,2/9 


3/8, 1/4, 5/24 










2 


(e) 


1/3,1/4,8/9 


-1/8,1/4, 7/8 










2 


(e) 


1/3, 1/4, 1/18 


7/8, 1/4, 1/244 










2 


(e) 


1/3,1/4, 7/18 


3/8, 1/4, 3/8 










2 


(e) 


1/3,1/4,5/9 


3/8,1/4, 13/24 


TAO-5 


6 Pm 


a 


(-1,0,0), (0,0, 2), (0,3,0) 


1 


(a) 


2/3,0,2/9 


5/6, 0, 1/4 


(12-atom) 






(0,0, 1), (-2, 0,0), (-1,-2,-1) 


1 


(a) 


2/3,0, 5/9 


0, 0, 7/12 










1 


(a) 


2/3,0,8/9 


2/3,0,11/12 










1 


{b) 


2/3,1/2,2/9 


5/6,1/2,1/4 










1 


(b) 


2/3,1/2,5/9 


1,1/2,7/12 










1 


(b) 


2/3,1/2,8/9 


2/3,1/2,11/12 










2 


(c) 


1/3,1/4,1/9 


1/2,1/4,1/12 










2 


(c) 


1/3,1/4,4/9 


1/6, 1/4, 5/12 










2 


(c) 


1/3, 1/4,7/9 


1/3, 1/4,3/4 


TAO-6 


11 P2i/m 


a 


(1,0,0), (0,0,1), (0,-6,0) 


2 


(e) 


1/3, 1/4, 13/18 23/24, 1/4, 17/24 


(12-atom) 






(0,0,-1), (-1,0,0), (2,4, -5) 


2 


(e) 


1/3,1/4,2/9 


11/24, 1/4,5/24 










2 


(e) 


1/3, 1/4,8/9 


5/8, 1/4, 7/8 










2 


(e) 


1/3, 1/4, 1/18 


7/24, 1/4, 1/24 










2 


(e) 


1/3,1/4, 7/18 


1/8,1/4, 3/84 










2 


(e) 


1/3,1/4,5/9 


19/24, 1/4, 13/24 



TABLE II: Pathways for the a— >a; phase transition in Ti. Each pathway is defined by a maximal common subgroup of both 
the a and u) phases and Wyckoff positions. A supcrccU containing both phases is given in terms of the lattice vectors of a 
and uj unit cells; the lattice vectors determine [n] and [m]. The atomic positions for each phase are given in the setting of the 
common subgroup; together the positions define Uj and 5j. The multiplicity m for each Wyckoff position is given with the 
label of the position. 
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Silcock TAO-1 2 3 4 5 _6_ 
Homogeneous barriers (in meV/atom) 
Elastic: 3.7 18 21 3.7 3.7 21 74 

TB Landscape: 60 43 61 83 80 84 81 
TB NEB: 54 24 52 - - - 

Ab initio NEB: 31 9 58 32 68 37 69 

Transformation information 
Supcrccll size: 12 6 6 12 12 12 12 
Orientation II I II II II II II 

Relations: I: (OOOl)a || (0111)^,, [1120]a || [0111]^ 
II: (0001)„ II (1120)^, [1120]c || [OOOlj^ 



TABLE III: Comparison of lowest landscape barrier path- 
ways. Energy barriers: Four different methods for calculating 
the energy barrier for the three pathways are shown, from 
least accurate to most accurate. The elastic barrier only ac- 
counts for the strain in each pathway. The landscape barrier 
uses a simple combined shuffle for each, and a tight-binding 
total energy. Finally, the NEB calculation is done with the 
tight-binding method for the lowest three and ab initio for all 
seven to accurately determine the barrier. Orientation Rela- 
tions: The relative orientation between a and co is shown for 
each pathway. N.B.: Silcock and the TAO-2 through TAO-6 
pathways have the same orientation relations. 




Ij 



o 




(OOOl)-a plane 



(Olll)-a; plane 



FIG. 2: TAO-1 acting on a single a basal plane. A hexagon in 
the basal plane is broken into two 3-atom pieces. Each 3-atom 
piece swings out of the basal plane in opposite directions, to 

create half of the honeycomb lattice in two parallel (0001) 
planes of uj. The remaining atom in the center forms the A 
sublattice lying in the (0001) planes of cu. The next a basal 
layer creates the dashed atoms in the lo honeycomb. This 
transformation turns the (0001) plane of a into the (0111) 
plane of u>. 



TAO-1 has the lowest energy barrier that best trades 
off between shuffle and strain. During the transforma- 
tion, the closest nearest neighbor distance is 2.63 A, 
which is larger than the 2.55 A value for TAO-2, 2.58 A 
for Silcock, TAO-3, TAO-4, and TAO-6. Only TAO-5 
has a comparable closest nearest neighbor distance of 
2.63 A; this larger distance helps to explain the lower bar- 
rier compared to TAO-2 with the same strains. Because 
the nearest-neighbor distances in a and lu are 2.93 A and 
2.65 A, respectively, it is not surprising that TAO-1 has 
the lowest barrier. The larger barrier of TAO-5 compared 
to TAO-1 may be due to more complicated geometric 
changes during the pathway. The lack of a simple, sin- 
gle criterion to explain the difference in energy barriers 
shows the necessity of ab initio studies. 



VI. CONCLUSION 

We have presented a general systematic pathway gen- 
eration and sorting algorithm for martensitic transfor- 
mations. The method enumerates all possible pathways 
within a few geometric restrictions. When applied to the 
Ti a^u transformation, a previously unknown pathway 
emerges with a barrier much lower than all other path- 
ways. Geometric criteria are useful in reducing the large 
set of possible pathways to a more manageable set, but 
ab initio studies arc required to ultimately find the lowest 
energy barrier pathway. 
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